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Abstract 

We study the properties of an anisotropically paired superconductor in the presence of a specularly 
reflecting surface. The bulk stable phase of the superconducting order parameter is taken to have 
d^2_y2 symmetry. Contributions by order parameter components of different symmetries vanish in 
the bulk, but may enter in the vicinity of a wall. We calculate the self-consistent order parameter 
and surface free energy within the quasiclassical formulation of superconductivity. We discuss, in 
particular, the dependence of these quantities on the degree of order parameter mixing and the 
surface to lattice orientation. Knowledge of the thermodynamically stable order parameter near a 
surface is a necessary precondition for calculating measurable surface properties which we present 
in a companion paper. 

To appear in J. Low Temp. Phys., Vol. 101, Dec, 1995 



I. INTRODUCTION 

Anisotropic superconductors react sensitively to electron scattering from impurities, interfaces and surfaces [0. 
Scattering leads, in these cases, to pair-breaking and a modification of the superconducting state. This paper is a 
first step in a systematic study of the effects of surface scattering on the superconducting properties of a layered 
system with d-wave pairing. The possibility of d-wave pairing is currently under serious consideration as a promising 
candidate for the proper order parameter symmetry in the high-Tc cuprates |^-^ . The presence of a d-wave pairing 
interaction results in a strongly anisotropic superconducting order parameter, and one expects that surfaces will exert 
a sizable infiuence on the pair condensate near the boundary p|,P-pO[| . Correspondingly, the physical properties of a 
d-wave system near an interface may be profoundly different from those of an s-wave superconductor. The distinction 
is primarily due to the strongly enhanced pair-breaking effects in a d-wave system. Many standard experiments, 
such as tunneling measurements or electromagnetic absorption at surfaces (in the very long mean free path limit), 
are particularly sensitive to the superconducting state near the interface. This suggests that experiments directed at 
measuring effects deriving from the interaction of the superconducting order parameter with a surface or interface 
may serve as reliable probes for the symmetry of the superconducting order parameter. The optimal surfaces for 
studying d-wave symmetry are those oriented perpendicular to the layers, and we consider this geometry exclusively 
here. 

Because the superconducting order parameter is itself not directly accessible to experiments, any evidence for an 
unconventional symmetry of the order parameter must be more or less indirect. Theory is required to extract infor- 
mation on the order parameter from the experimental properties of the superconductor, and to point out measurable 
effects which reflect clearly the symmetry of the underlying order parameter. This paper and the companion paper, 
[II] |2l|, discuss specific aspects of anisotropic superconductivity near surfaces. We consider systems with tetrag- 
onal (i^4fi) crystal symmetry, which is a good approximate symmetry of the cuprate superconductors, and focus 
on a non-degenerate order parameter (Cooper pair wavefunction) of even parity. A group theoretical classification 
of superconducting order p2| leads, in this case, to four different types of bulk order parameters corresponding to 
the four one-dimensional, even parity representations of the D4/j group: the Ai, or identity representation, the A2 
representation, and the Bi and B2 representations. Whichever representation leads to the largest attractive cou- 
pling constant for Cooper pairs will determine the symmetry and transition temperature of the bulk superconducting 
state. A surface, on the other hand, may break the D^ih symmetry and thus may admit local contributions from all 
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other representations, to a degree depending on the magnitude of their coupHng constants. In addition, surface pair- 
breaking may suppress the dominant bulk order parameter and generate surface states of a new symmetry. Surface 
pair-breaking and the mixing in of new symmetries thus leads to measurable modifications of the superconducting 
properties near a surface, which carry information on the symmetry of the superconducting order parameter. 

We investigate both the thermodynamic and spectral properties of a superconductor in the vicinity of a specular 
interface. We first focus on calculations of the spatial dependence of the order parameter (gap function) and of the 
surface free energy as functions of the relative surface to crystal lattice orientation. We take the bulk order parameter 
to be of Bi d-wave symmetry as predicted by microscopic models [^|J^], and indicated by various experiments. In 
the presence of a surface, we consider three main order parameter categories: those possessing a single Bi d-wave 
component (all other coupling constants being assumed negligible), those possessing a linear combination of Bi and 
B2 d-wave components (the Ai and A2 coupling constants being assumed negligible), and those possessing a linear 
combination of the Ai, A2, Bi, and B2 components (i.e. all components allowed by symmetry may couple in). We 
discuss to what extent the mixing of these various symmetry components affects the overall order parameter structure 
in the vicinity of a surface as well as the corresponding effect on the surface free energy. We use these results in [II] 
to study the influence of surface pair-breaking on the tunneling density of states. 

We perform all of our calculations within the quasiclassical formulation of superconductivity. This is a fully 
self-consistent theory which requires the numerical calculation of the quasiclassical propagator from which the ther- 
modynamic and spectral data may then be extracted. The formulation is justified as long as the superconducting 
coherence length is large compared to the Fermi wavelength. This approximation is known to be excellent for tradi- 
tional superconductors, but it may be of reduced validity for the high-Tc cuprates due to their substantially shorter 
coherence lengths. In section II we give a brief review of the quasiclassical framework, and in section III we present 
our numerical results for the gap-function and free energy. A discussion and summary are contained in section IV. 



The quasiclassical formulation of superconductivity was first stated by Eilenberger ||2^, Larkin and Ovchinikov 
|p^ , p5| , and Eliashberg |^ . This theory is capable of describing both the equilibrium and dynamical properties of 
conduction electrons in either the normal or superconducting state and may be understood as a generalization of 
Landau's theory of normal Fermi liquids. In this section we give a brief summary of only those equations relevant to 
the current analysis. Readers interested in a more detailed presentation are referred to the original papers mentioned 
above, and to the recent review articles p^-|30[|. As our primary resource regarding the quasiclassical formulation we 
turn to reference |^8|] , whose notation we adopt whenever possible. 

In this paper we present calculations of the thermodynamic properties of a superconducting system in equilibrium 
and, correspondingly, require only the quasiclassical theory in its imaginary-energy (Matsubara) representation. The 
fundamental quantity in the quasiclassical formulation is the matrix propagator for quasiparticle excitations, g. For 
spin-independent systems the propagator has the simple form: 



which is a 2 X 2 matrix in particle-hole (Nambu) space. The propagator ^ is a function of a momentum on the 
Fermi surface, P/, a real-space position, R, and an imaginary excitation energy, ze„. The diagonal components of g 
carry information on the spectrum of Bogoliubov quasiparticles, while the off-diagonal components are related to the 
Cooper pair amplitude. The functions / and g are related to / and g by fundamental symmetry identities [ p8| . 
The central equation obeyed by the propagator g is the quasiclassical transport equation 



II. QUASICLASSICAL THEORY 




(1) 



ienh - A(p/,R) , 5(p/,R;e„) (p/) • VR5(p/,R;e„) = 0, 



(2) 



supplemented by a normalization condition 




^Note that we are using the terms pairing self-energy, order parameter, and gap function interchangeably. 
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A (p;, R) = / dp} n{p'f )V (p/, p}.) / (p'^, R; e„) 



(4) 



Here V |^p/,Pjj is the pairing interaction, which determines both Tc and the symmetry of the order parameter, 

and /(p/,R;e„) is the off-diagonal part of the full matrix propagator 5(pj^,R;e„). The integration over py is to 
be interpreted as a normalized Fermi surface integral such that § dpf 7i(p/) — 1, where ?^(p/) is the Fermi surface 
anisotropy factor (see [II]). Equations (^) and (^) must be solved self-consistently for the position dependent gap 
function A (p j , R) . 
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Fig. 1. This figure defines the geometrical configuration and the central descriptive parameters of our calculation. 
We sketch the clover-shaped Bi d-w&ve gap structure whose bulk orientation is fixed to the crystal lattice. The 
angle (j)o is the surface to lattice orientation angle {(f)o = 45 deg in this case), while ^out defines the trajectory of the 
out-going quasiparticle. The trajectory of the incoming quasiparticle is specified by {(j>^^ — (j>o) = tt — ((/>out ~ 't'o)- 



We consider a semi-infinite system occupying the space x > with a single specularly reflecting interface located 
at a; = 0; this configuration is depicted in Fig. 1. The boundary condition for the quasiclassical propagator ^ at a 
specular wall requires continuity along a reflected trajectory: 

5(P/,in'^^waU'^") = 5(P/,out' ^^waU' ^"); (5) 

where the out-going momentum vector, Pyout' given in terms of the in-coming momentum vector, Py jj^, and the 
surface normal, fi, by 



P/,out ~ P/,in 




Note that an anisotropic order parameter need not, in general, be continuous at the surface along a given trajectory 
since the momentum vector abruptly changes direction there. 

In our model we assume a homogeneous order parameter in the plane of the surface; in this way there is only one 
relevant spatial degree of freedom which allows us, in the following, to replace R by x. Further, we consider layered 
superconductors with negligibly small interlayer coupling. We model these systems by a cylindrical Fermi surface with 
no dependence of the gap function on pz so that the Fermi surface vector p f can be replaced simply by the azimuthal 
angle 0. These simplifications allow us to rewrite the quasiclassical transport equation in the form 

ie„f3 - A{(p,x) , g{(j),x;en) 

which amounts to an ordinary first order differential equation along a given "trajectory" (incoming and outgoing) 
which is defined by the parameter (p. These trajectories may be understood as "classical" trajectories of particles 
moving with velocity Vf along the direction </> and being reflected at the surface. A typical trajectory is shown in 
Fig. 1 where we have chosen to measure the angle (p relative to the crystal a-axis. In our simplified model the Fermi 
surface parameters are assumed isotropic in the a6-plane, and the reduction to tetragonal symmetry comes about only 



-I- ivf cos 4>-^gi<j), x; e„) = , 



(7) 
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via the pairing interaction. The use of more general Fermi surface parameters would not alter the main results of our 
calculations.^ 

We incorporate contributions to the pairing self-energy from the Ai, A2, Bi, and B2 representations of the tetragonal 
group. Our approximate order parameter is taken to be a linear combination of simple trigonometric functions of the 
angle 

= As{x) X 1 + Agi (x) X cos 4(/) + Ag2 (x) x sin ic/) 

+Adi (x) X cos 24) + Ad2 (x) X sin 20 (8) 

In the usual notation, 1 is an s-wave function, cos4</) and sin4(/) are g-wave functions of Ai and A2 symmetry 
respectively, and cos2(/) and sin 20 are d-wave functions of Bi and B2 symmetry. In accordance with equation (^), we 
decompose the gap equation into the following five self-consistency equations: 

A,(a;) =y.r|]^^/(0',:r;e„), (9) 

Ag, (</.,x) = V^g,r|]^^2cos40cos407(0',a;;e„), (10) 

Ag,{cl,,x) =14,r|]^^2sin4</)sin4</)7(</)',a;;e„), (11) 

Adi(</>,x) =l/d,r|]^^2cos20cos20'/(0',x;e„), (12) 

Ad, (0,x) =Fd,r|]^^2sin20sin20'/(0',x;e„). (13) 

where the Vx are coupling constants, and the indices s, d, and g refer to s-wave, d-wave, and gr-wave respectively. 
Note that As{4>,x) and Ag-^{(j>,x) have Ai symmetry, Ag2(0, x) has A2 symmetry, and Adj^(0, x) and Ad2{(f',x) have 
Bi and B2 symmetry respectively. One can always eliminate the coupling constants Vx , along with the frequency-sum 
cutoffs eco, in favor of the transition temperatures via the customary BCS relation Tcx ~ l-lSccoe"^/^-^ , where X 
denotes the subscripts s, gi, 52, ^i, and d2- In general, Vx may be positive (attractive) or negative (repulsive). 

Equations (0)-(^ must be solved self-consistently for the spatially dependent gap amplitudes, Ax{x), which are 
the upper right elements of the corresponding pairing self-energies. 

We may then proceed to calculate various thermodynamic quantities such as the free energy difference per unit surface 
area: 

/•OO 

^surf^ / dx {T{x) - J^tuik} , (15) 
Jo 

where J-{x) is the free energy density in the presence of the wall, 



T{x)=NiEF){ln[^] \AAx)\' + ^ 



ln(-^) \Ag,{x)\^+lnf^]\Ag.M\ 

^cgi / V J cg2 , 



significant influence of the Fermi surface geometry is expected if the Fermi energy lies near a Van Hove singularity (for a 
discussion of these topics see Refs. p]]-|3^). The effect of non-trivial Fermi surfaces will be discussed in a future publication. 

^The precise ^-dependence of the order parameter is still controversial [^. For this study we decided to parameterize the </) 
dependence in a straight forward way. An alternative approach would have been to take the order parameter from calculations 
in the single-exchange graph approximation. M. 



4 



In 



T 



cdi 



|Ad,(x)|2+ln 



/ T 



\Tcd2 



+ iN{EF)k^T J2 (tr4 [f3 (g ( 



,x;X) -go((/>;A))] -2i 



.^|A(0,x)|2 



A2 



(16) 



and Tbuik is obtained by setting A((/), a;) = Abtiifc((/)). Here go is the normal state propagator, and A(0, x) is the fuU, 
(^-dependent gap ampHtude. In the next section we make use of the preceding equations to calculate the various gap 
amplitudes and the surface free energy for various surface orientations. 



III. RESULTS 



We consider a superconductor-insulator interface, with the superconductor occupying the space a; > 0. We compute 
the constituent quantities as a function of x out to a distance of ~ 30^ = Hvf /TrAmaxix = oo)) from the wall. We 
consider an ideal interface by imposing a specular boundary condition at the wall, which amounts to the condition: 

9 (0in' a; = 0; e„) = g {(pont^x = 0; e„) , (17) 

where (/>out complementary angle to (obtained upon specular reflection of the trajectory) and is given by 

('/'out - 0o) = TT - {(f)^^^ - (t)o). 

We divide the presentation of our results into two subsections. In the first subsection we consider a simple d-wave 
order parameter of the form A((j),x) = A^^^ (a:) cos 2(/), and in the second we include the possibility of admixtures of 
components of different symmetries which are induced by the interaction with the surface. 



A. Simple d-wave Model 



In this section we consider a pairing interaction of the form V{(l),(j)') — 2Vdi cos2(/)COs2(/)', which corresponds to 
a purely d-wave order parameter of Bi symmetry. This is the symmetry favored by microscopic theories to give 
the highest Tc ||,^. The order parameter has the form: A((/), a;) = A^j (a;) cos 20 where Ad;^{x) must be evaluated 
self-consistently from equations (|^) and (^2|) . In Fig. 2 we plot the gap amplitude A^i (x) as a function of x for various 
surface to lattice orientation angles 0o (the angle 0o is defined in Fig. 1). Note that for an orientation angle (po — the 
gap amplitude is not suppressed by the presence of the wall. This is analogous to the isotropic s-wave case. However 
as the surface to lattice orientation angle increases the gap amplitude rapidly becomes suppressed in the vicinity of 
the wall and it completely vanishes at the wall for 0o ~ 45 deg. Note that this suppression heals on the scale of about 
six coherence lengths. 
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Fig. 2 The order parameter amplitude as a function of position when only one d-wave component is present, for 
various surface to lattice orientation angles, (\)o- 
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Fig. 3 The local order parameter experienced by a quasiparticle moving along a trajectory with angle (j^out f^'" * 
surface to lattice orientation angle of 4>o = 30 deg. The parameter p is the coordinate position along the trajectory 
and is related to x via p = x/ cos 0out ~ 4'o- The point of reflection is taken to be p = 0. 
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Fig. 4 The surface free energy as a function of surface to lattice orientation, 0o, for various degrees of mixing in 
units of (A''(_Ef)| Aj.j^j^^_jj^j]^(r = 0)|^^) . The two center lines correspond to the low transition temperature regime 
{Tcx/Tc = 0.2) while the lower two linos correspond to the high transition temperature regime (Tcx/Tc = 0.6). 
The solid line corresponds to purely Bi pairing. 



The angular dependence of the total gap suppression can be understood within the framework of the quasiclassical 
approximation. A quasiparticle moving along a classical trajectory experiences the local order parameter with mo- 
mentum along (j) (i.e. in the direction of the trajectory). For a surface parallel to one of the crystal axes, the local order 
parameter is constant along the trajectory because A^^ (^i^) — ('/'out) ^^'^ solution of the transport equation 
along this trajectory is identical to the corresponding solution in the bulk. As a consequence, the self-consistently 
determined order parameter is identical to the bulk order parameter for this orientation of the surface and the surface 
free energy is zero. For any other surface orientation one has Adj((^j^) 7^ Adi((/'out)' propagator /, and 

consequently the self-consistent order parameter, will be modified by the surface. In Fig. 3 we plot the local order 
parameter A as a function of the spatial position, p, along the trajectory for a surface orientation = 30 deg and 
several different trajectories. The point of reflection is taken to be at p = 0, with negative p values denoting an 
incoming particle. This depicts the local gap amplitude that a quasiparticle experiences along its trajectory. As 
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previously mentioned, a discontinuous jump in the order parameter upon reflection at the wall can be observed along 
all trajectories (except when (/lout ~ '^o = deg) due to the discontinuous change in momentum direction. Parenthet- 
ically, we note that an asymptotic change in sign of the order parameter along a trajectory (e.g. Fig. 3, trajectories 
'^out ~ 0o = 30 deg and 45 deg) can give rise to significant spectral effects such as zero-energy bound states |]37| and 
Tomasch oscillations pj] ]. Our discussion of spectral effects will appear in a companion paper. 

The surface free energy for the case of a simple d-wave order parameter is displayed as a function of the surface 
orientation angle (po in Fig. 4 as the solid line. The surface energy is zero, as expected, for a surface parallel to one 
of the crystal axes (^o = Odeg, 90 deg, . . . ). For any other orientation the surface energy is positive as a consequence 
of the surface pair-breaking. Observe, in this case, that the energy values are very strongly dependent on the relative 
orientation of the gap nodes with respect to the wall. These qualitative features of the pure d-wave model will persist 
even if components of other symmetries are mixed in. However, as we discuss in the following subsection, mixing in 
leads to significant quantitative differences in the order parameter structure and the surface free energy. 




Fig. 5 The individual gap amplitudes, Ax{x), as functions of x for the low coupling constant regime {Tcx /Tc = 0.2) 
at a surface to lattice orientation angle of (j)o — 30 deg. 




Fig. 6 The individual gap amplitudes, Ax(x), as functions of x for the low coupling constant regime {Tcx/Tc = 0.6) 
at a surface to lattice orientation angle of (po ~ 30 deg. 
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B. Extended d-wave Model 



In this section we consider an admixture of the previously mentioned components Ag, A^^, Ag^, and Ag^ to the 
dominant A^^ component. Thus the total gap amplitude has the form 

A(^, x) = As{x) + Agj (x) cos 40 + Ag^ (x) sin 40 

+ Adj (a;) cos 20 + Ad2 (a;) sin 20. (18) 

The degree to which a component Ax mixes in is determined by the corresponding transition temperature, Tcx- The 
dominant component has the highest transition temperature {Tcdi ) and provides the scale for all other temperatures 
in the calculation. Accordingly, wc have Tc = Tcdi and then the ratios Tcx/Tc become the parameters of the 
calculation. We consider only coupling constant and temperature values for which the non-dominant components 
(i.e. Tf-x/Tf, < 1) vanish in the bulk, but may appear near the wall. The admixture of the subordinate components 
changes the momentum dependence of the gap function. For instance, an admixture of A^^ causes a rotation of the 
total gap (e.g. Ad^ alone is equivalent to a 45deg rotation in momentum space). The effective gap maximum and 
the angular position of the gap's nodes thus vary, in general, with distance from the wall. By means of mixing in the 
new components, the system is able to alter the gap's angular and spatial structure and thereby further lower its free 
energy. 
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Fig. 7 The maximum of the total order parameter as a function of x in the low coupling constant regime, 
Tcx/Tc = 0.2, for various surface to lattice orientation angles, <f>o- 
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Fig. 8 The maximum of the total order parameter as a function of x in the low coupling constant regime, 
Tcx/Tc = 0.6, for various surface to lattice orientation angles, 0o. 
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The individual gap amplitudes arc plotted versus distance from the wall in Figs. 5 and 6 for a temperature T/Tc = 
0.7. These two graphs portray the behavior of the gap amplitudes for two different coupling constant regimes. For 
simplicity we have set all of the admixture components to have the same value for Tcx jTc- Note that in both cases the 
two d-wave components enter with by far the largest amplitudes and are truly the dominant representatives. Indeed, 
a computation without the s and g-wa.ve components left the amplitudes of the two d-wave components virtually 
unaltered. 

The corresponding maximum gap values are plotted versus distance from the wall in Figs. 7 and 8 and should be 
compared with Fig. 2 where the admixture is absent. Clearly the system is now able to compensate and avoid (to 
some degree) gap reduction with the level of compensation being greater for larger Tcx ■ Note that for both transition 
temperature regimes the total gap maximum for a surface orientation of = 45 deg is nearly unaltered from the 
corresponding curve without mixing. This is a special symmetric case whose physics depends strongly on the ratio 
T/Tcx and thus presents a potential mechanism for actually measuring the admixture coupling constants. 
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Fig. 9 The angular orientation of the order parameter maximum with respect to the surface normal as a function 
of X for three different degrees of mixing: a single Si d-wave component, a Bi and a B2 d-wave component, and 
full mixing. 
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Fig. 10 The maximum of the total order parameter as a function of x for three different degrees of mixing: a 
single Bi d-wave component, a Bi and a B2 d-wave component, and full mixing. 
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Another striking quality of the self-consistent solution is the complete reorientation of the total gap's nodes toward 
the optimal configuration (i.e. A((/)jj^) ~ A((/)q^^)) at the wall, independent of the surface orientation and even for 
surprisingly small admixture coupling constants. This effect can be plainly seen in Fig. 9 which shows the angular 
orientation of the maximum of the gap function as a function of x for different degrees of mixing. In the bulk 
(pmax = Odeg, ±90deg,. . . , however in the vicinity of the wall the total gap rotates in an attempt to orient the 
maximum gap along the surface normal {(^^max — 4>o = deg, ±90 deg,. . . ) at a; = 0. This overall rotation of the 
bulk Bi d-w&ve gap function is achieved by the mixing in of the B2 d-wave component. The mixing in of the other 
components can affect the magnitude and the shape of the total gap function, but they may not bring about an 
overall rotation. This is presumably the reason for the system's preference of the B2 d-wave component over the other 
symmetry components. 

The maximum gap values for a given surface orientation are plotted as a function of x in Fig. 10 for various degrees 
of mixing. Observe that in the immediate vicinity of the wall the addition of the second (-B2) d-wave component 
has a substantial effect, while the addition of the remaining three has a less pronounced effect. This interpretation 
is substantiated by the results of the free energy calculations which are plotted in Fig. 4. We display results for a 
single d-wave component, two d-wave components, and the fully mixed case. Clearly, the inclusion of the Ai and A2 
components in addition to the Bi and B2 components improves the free energy very little. Indeed, the lines lie within 
a few percent of each other for most surface orientations. 
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Fig. 11 The orientation and magnitude of the total order parameter as a function of distance from the wall 
for a surface to lattice orientation angle of 4>o = 30 deg. The bulk order parameter (at x = 00) is included for 
comparison. 
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Fig. 12 The local order parameter experienced by a quasipaxticle moving along a trajectory of angle 0out * 
surface to lattice orientation angle of <po = 30 deg. This is the case of an order parameter consisting of Bi and 
B2 d-wave components. The parameter p is the coordinate position along the trajectory and is related to x via 
p = x/ cos 0out ~ 'I'o- The point of reflection is taken to be p = 0. 
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The crucial importance of the B2 components is emphasized in Fig. 11, which depicts the angular orientation of 
the gap function with respect to the surface normal as a function of position. In this figure we consider an order 
parameter consisting of only the Bi and B2 d-wavc components. The addition of B2 allows the system to rotate the 
total gap maximum toward its optimal orientation (normal to the surface) at the wall. This occurs for both large 
and small coupling constants and represents one of the most prominent results of this calculation. We fully expect 
that this could have profound consequences for measurements since even a small admixture assures that the order 
parameter assumes a symmetric orientation at the wall regardless of its bulk orientation. 

We conclude our presentation of the numerical results by presenting a plot of the local gap experienced by a 
quasiparticle along a given trajectory. The results for the Bi + B2 mix arc shown in Fig. 12 and should be compared 
with Fig. 3. The fact that the total gap has rotated toward its optimal orientation at the wall can be seen, for example, 
by examining the <?i>out ~ "Ao = trajectory where the local gap is now greater at the wall than in the case without 
mixing (compare with Fig. 3). 

IV. CONCLUSIONS 

In this paper we have examined the thermodynamic properties of an anisotropically paired superconductor in the 
vicinity of a perfectly reflecting surface. In particular we considered an order parameter which can be expressed as 
a linear combination of components transforming as the Ai, A2, Bi, and B2 representations of the I?4^ (tetragonal) 
group. We carried out self-consistent calculations of the various gap amplitudes within the quasiclassical theory of 
superconductivity. Such calculations are needed as a starting point for calculations of a number of physically relevant 
quantities such as the excitation spectrum, the surface impedance, Josephson couplings, etc. Another important 
example is the tunneling density of states which is discussed in [II]. 

Starting from the assumption that a d-wave order parameter transforming as the Bx representation is the bulk 
stable phase, we showed that even for a small T^b^ the mixing in of a B2 component can have a profound effect on the 
order parameter structure near the wall. This arises from the consequent effective rotation of the total gap toward the 
optimal (i.e. symmetric) orientation in the immediate vicinity of the wall. The effects of the mixing in of the Ai and A2 
components has a smaller, yet noticeable, influence on the gap function, while their effect on the surface free energy is 
considerably less. Thus the presence of the Ai and A2 symmetries has a minor effect on the thermodynamic properties 
of the system, but can have a more pronounced effect on the excitation spectrum. The strongest pair-breaking for a 
Bi symmetry order parameter is obtained at a 45 deg surface, and for a B2 order parameter at a surface along one of 
the crystal axis (Odeg surface). We considered here ideal surfaces without roughness or degradation. This represents 
the optimal configuration for probing the order parameter by surface measurements. Surface imperfections will wash 
out, to some degree, the measurable orientational effects of an anisotropic order parameter. As a resiflt, experiments 
should be performed preferentially at clean and well characterized surfaces oriented perpendicular to the a6-planes. 

Note Added in Proof. After the submission of this manuscript we received a preprint by M. Matsumoto and 
H. Shiba in which the possibility of the mixing in of subdominant symmetry components which break time-reversal 
symmetry is discussed. We have not allowed for the possibility of such states in our calculations. These time-reversal 
symmetry breaking states are stable in certain regions of the coupling-constant /temperature parameter space. All of 
the calculations in the current manuscript, however, were done in a region of parameter space where these time- reversal 
breaking states are not stable. 
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